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Summary 

This  document  describes  the  fundamentals  of  heat 
transfer  measurement  techniques  in  high  enthalpy 
flows.  After  description  of  basic  relations  for  the 
stagnation  point  heat  flux  rate  to  a spherical  surface 
in  a high  enthalpy  flow  field,  different  heat  flux 
measurement  techniques  are  discussed.  Attention  is 
paid  to  design  aspects,  mathematical  model  for  the 
heat  flux  rate  evaluation  and  application  limits  of 
each  sensor  type,  in  order  to  provide  a useful  and 
practical  script  for  the  reader  with  respect  to  the 
choice  of  an  adequate  sensor  type  for  different 
requirements.  Sensor  calibration  and  comparative 
measurements  using  different  heat  flux  sensors  in 
the  arc  heated  facility  LBK  are  described  in  the  last 
two  chapters. 

1.  List  of  symbols 

A cross  section  area 

C specific  heat 

Dn  binary  diffusion  coefficient  for  atoms 

and  molecules 
h enthalpy 

hD  enthalpy  of  dissociation  per  unit 

mass  for  the  gas  in  the  external  flow 
K thermal  conductivity 

LBK  arc  heated  facility  of  DLR 

Le  Lewis  number 

m mass  flow  rate 


Nu  Nusselt  number 

P pressure 

Pr  Prandtl  number 

q heat  flux  rate 

Re  Reynolds  number 

s wall  thickness 

T temperature 

u flow  velocity  component  in  the  direc- 

tion parallel  to  the  body  surface 

£ emissivity 

P density 

Subscripts 

c 
e 
s 
w 
0 
1 
2 
00 

2.  Introduction 

Intrusive  measurement  techniques  in  supersonic 
and  hypersonic  flows  have  been  a challenging  task 
for  engineers  for  decades.  Recent  progresses  in 
non-intrusive  measurement  techniques  and  numeri- 
cal codes  have  led  to  a better  understanding  of  local 
flow  phenomena  and  to  further  improvement  of 
intrusive  measurement  techniques  in  aerodynamics. 
Because  of  its  complexity  the  measurement  of  the 


coolant 

boundary  layer  edge 
stagnation  point 
wall 

stagnation  condition 
inlet  side 
outlet  side 
ambient,  free  stream 
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heat  flux  rate  remains  one  of  the  key  techniques, 
which  have  to  be  improved  with  respect  to  meas- 
urement accuracy  and  repeatability.  The  heat  flux 
rate  is  usually  determined  from  the  measured  tem- 
perature development  on  the  surface  or  inside  the 
probe  and  applying  a mathematical  model  with 
some  assumptions.  Any  deviation  from  a perfect 
contact  between  the  temperature  sensor  and  probe 
or  other  satisfaction  of  assumptions  leads  to  inaccu- 
racies. 

Heat  flux  determination  in  high  enthalpy  flow  fields 
is  accompanied  with  some  additional  effects.  Since 
the  flow  field  of  most  of  the  high  enthalpy  facilities 
is  in  thermo-chemical  non-equilibrium,  the  gas  sur- 
face interaction  phenomena  influence  the  heat  flux 
rate  significantly.  Surface  catalysis  and  emittance, 
which  are  a function  of  the  surface  temperature, 
play  an  important  role  in  such  environments.  The 
boundary  flow  around  the  probe  is  also  influenced 
by  these  effects  and  there  is  a continuous  coupling 
between  the  flow  and  model  structure.  Differences 
in  the  operation  of  different  high  enthalpy  facilities 
give  rise  to  different  requirements  on  the  heat  flux 
sensors.  Short  duration  high  enthalpy  facilities  like 
shock  tunnels,  impulse  facilities,  etc.  have  a very 
short  testing  time  of  several  milliseconds.  Therefore 
only  heat  flux  sensors  with  a very  short  response 
time  can  be  used  in  these  facilities.  Because  of 
short  flow  exposure  time  of  the  sensor  the  surface 
temperature  is  low.  Arc  heated  and  induction 
heated  facilities  have  a longer  testing  time,  which 
would  lead  to  higher  surface  temperature  of  the 
stationary  and  uncooled  probes.  The  practical  reali- 
sation of  the  cooling  and  insulation  of  the  probe 
influences  the  measurement.  Another  option  is  the 
application  of  uncooled  transient  heat  flux  probes, 


whose  reliability  depend  on  some  further  parame- 
ters like  sweep  speed,  flow  homogeneity,  etc.. 

Because  of  above  mentioned  diversity  in  the  re- 
quirement on  heat  flux  sensors  different  types  of 
measurement  techniques  will  be  described  and 
compared  in  this  document.  Since  some  techniques 
can  be  used  in  both  short  and  long  test  duration 
facilities,  the  layout  of  this  script  is  chosen  in  such 
a way  that  similar  techniques  are  described  in  the 
same  chapter. 

3.  Heat  flux  at  high  surface  temperatures 

The  heat  transfer  process  at  hypersonic  velocities  is 
complicated  by  high  temperature  effects,  which  are 
not  present  at  low  speed  flights.  The  shock  wave  in 
front  of  the  critical  components  of  a hypersonic 
vehicle  like  nose  cap,  wing  leading  edge,  etc.  heats 
air  to  high  temperatures  at  which  dissociation  and 
ionization  of  the  gas  species  take  place  (Fig.  1). 


Figure  1.  Stagnation  region  flow  field 

Since  the  dissociation  and  ionization  processes 
proceed  at  finite  rates  and  the  relaxation  length 
behind  a bow  shock  is  small  the  thermochemical 
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equilibrium  is  not  necessarily  be  achieved  in  the 
stagnation  point  region.  Depending  on  the  shock 
strength,  i.e.  flight  Mach  number,  the  gas  at  the 
boundary  layer  edge  consists  of  different  composi- 
tion of  molecules,  atoms  and  ions.  The  vibrational 
rotational  and  translational  excitations  of  molecules 
are  not  in  equilibrium.  Due  to  viscous  forces  in  the 
boundary  layer  the  gas  is  decelerated  and  recombi- 
nation reactions  are  initiated.  Since  atoms  and  ions 
recombine  with  a high  specific  energy  release  the 
heat  transfer  rate  is  enhanced.  If  a fully  thermo- 
chemical equilibrium  of  the  gas  is  not  achieved  in 
the  boundary  layer,  the  heat  flux  rate  to  the  body 
depends  on  the  catalysis  of  the  surface. 


specific  heat 

Du: 

binary  diffusion  coefficient  for  atoms 
and  molecules 

h: 

enthalpy 

hD: 

enthalpy  of  dissociation  per  unit  mass 
for  the  gas  in  the  external  flow 

K: 

thermal  conductivity 

Qs- 

stagnation  point  heat  flux  rate 

u : 

flow  velocity  component  in  the  direc- 
tion parallel  to  the  body  surface 

x: 

radial  distance  along  the  body  surface 
from  the  stagnation  point 

P: 

viscosity 

p: 

density 

The  stagnation  point  heat  transfer  rate  to  a spherical 
body  in  a dissociated  viscous  boundary  flow  with 
thermochemical  equilibrium  was  described  by  Fay 
and  Riddell  [1,2]  as  follows: 


qs  =0.76  Pr~° 6 (pe[ie)° 4 (pwp.w)ai 


K ~hw) 


1 + (/.e0  52  - 1) 


rdue^ 
V dx  j 

fh  X 

11 D 

Kh%J 


•(1) 


The  non-dimensional  parameters  appearing  in  this 
formula  are  defined  as 

Prandtl  number 


Pr  = 


K 


Lewis  number 


(2) 


The  subscript  notation  is: 

0 : conditions  at  the  stagnation  point 

e : conditions  at  the  boundary  layer  edge 

w : conditions  at  the  body  surface 

The  stagnation  point  velocity  gradient  is  given  by 
Newtonian  theory  as: 

rdue\  =]_  [2 (pe-pj 
[dx)s  /?V  Pe 

where  R is  the  radius  of  the  sphere,  pe  and  px  are 
the  pressure  at  the  stagnation  point  of  the  boundary 
layer  edge  and  the  pressure  ahead  of  the  shock,  re- 
spectively. 

The  heat  flux  rate  of  a fully  catalytic  spherical  wall 
in  a frozen  boundary  layer  is  given  as 


Le  = 


P^-p^12 

K 


(3) 


The  other  notations  used  in  the  equation  (1)  are 
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qs  =0.76/V°6(pene)04  (p„n„,)' 
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Re  = 
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(8) 

(9) 


The  heat  flux  rate  to  a non-catalytic  surface  in  a 
frozen  boundary  layer  is  approximated  with  Le  = 0 
as: 


qi-0.76Pr-»-6(p,Mer(p^„)0' 


( due^ 

f i, 

1 _ _2_ 

dx  J 

S 

(6) 


The  equations  (1)  and  (5)  describing  the  heat  flux 
rate  to  an  arbitrary  surface  in  an  equilibrium  flow 
and  the  heat  flux  rate  to  a fully  catalytic  surface  in  a 
frozen  flow  field,  respectively,  are  essentially  the 
same.  They  differ  only  in  the  slightly  different  ex- 
ponent on  the  Lewis  number.  That  means  that  the 
surface  heat  flux  rate  is  nearly  the  same  whether  the 
flow  is  in  local  thermochemical  equilibrium  or  is 
frozen  with  a fully  catalytic  wall.  In  the  case  of  an 
equilibrium  flow  recombination  reaction  takes 
place  in  the  colder  region  of  the  boundary  layer 
itself,  releasing  chemical  energy,  which  is  trans- 
ported by  thermal  conductance  to  the  surface.  In  the 
frozen  boundary  layer  flow  recombination  occurs 
on  the  catalytic  wall  surface.  Figure  2 shows  the 

heat  transfer  coefficient  Nu/^[Re  as  a function  of 
the  recombination  rate  parameter  C,  [3].  Nusselt 
number  (Nu\  Reynolds  number  (Re)  and  C,  are 
defined  as: 


Nu  ■ 


<isCPj 


Kw(he-hw) 


(7) 


R is  the  universal  gas  constant. 

The  solid  lines  show  the  total  heat  flux  rate  for  both 
catalytic  (1)  and  non-catalytic  (2)  surfaces.  The 
heat  flux  rate  by  conduction  alone  to  a catalytic 
wall  is  shown  by  the  dotted  curve  (3).  Large  values 
of  C,  corresponds  to  an  equilibrium  boundary  layer 
flow.  At  this  range  the  total  heat  transfer  is  domi- 
nated by  conductive  heating.  As  mentioned  before 
in  an  equilibrium  flow  the  heat  flux  rate  to  both 
catalytic  and  non-catalytic  surfaces  is  the  same 
since  the  recombination  is  completed  in  the  bound- 
ary layer  flow.  In  contrast  the  heat  flux  rate  to  a 
non-catalytic  wall  depends  strongly  on  the  flow 
regime.  By  moving  from  high  C,  values  (equilib- 
rium) to  smaller  C,  values  the  flow  becomes  pro- 
gressively more  non-equilibrium  and  it  is  frozen  at 
the  left  end  of  the  curves.  It  can  be  seen  clearly  that 
the  heat  flux  rate  to  a non-catalytic  wall  (curve  2)  in 
the  frozen  flow  field  drops  by  a factor  of  more  than 
two  compared  to  the  equilibrium  flow.  The  differ- 
ence between  the  curves  1 and  3 represents  the  heat 
transfer  due  to  diffusion.  It  can  be  seen  clearly  that 
at  decreasing  Cl , i.e.  freezing  boundary  layer  che- 
mistry, the  conductive  heat  transfer  rate  to  the  wall 
is  decreasing.  Although  for  equilibrium  flows  the 
total  heat  transfer  is  essentially  all  conductive,  the 
diffusion  becomes  a larger  part  of  the  heat  transfer 
to  a catalytic  wall  in  strong  non-equilibrium  flow. 
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The  influence  of  the  non-equilibrium  phenomena  in 
the  gas-surface-interaction  in  high  enthalpy  flow 
field  has  to  be  considered  in  the  measurement  tech- 
niques of  the  heat  flux  rate  determination  in  high 
enthalpy  facilities. 


qg=qc+qr  (10) 

where  qg , qc  and  qr  represent  the  heat  flux  rate 

from  the  gas  to  the  surface,  conduction  inside  the 
model  and  thermal  radiation  of  the  model  surface, 
respectively.  The  thermal  radiance  is  a strong  func- 
tion of  the  surface  temperature 

ov  =s°(C-r„4)  <ii) 

where  s and  or  (=  5.67-1  O'8  [W/(m2  • K"4)])  are 
the  surface  emissivity  and  the  Stefan-Boltzmann 
constant,  respectively. 


Figure  2.  Heat  transfer  parameter  Ct  for  dif- 
ferent wall  catalysis  [1,2] 

4.  Heat  flux  measurement  techniques 

Fay  and  Riddell  equation  for  the  stagnation  point 
heat  flux  rate  does  not  include  radiation  effects. 
Since  the  radiance  of  a surface  is  a strong  function 
of  the  surface  temperature  and  its  development 
depends  also  on  the  thermal  properties  of  the  vehi- 
cle or  model  materials,  it  has  to  be  considered  by 
the  analysis  of  heat  balance  of  a solid  surface.  Al- 
though especially  in  high  density  gas  with  high 
temperature,  radiation  processes  may  affect  the  heat 
transfer  from  the  gas  to  the  model  significantly, 
thermal  radiation  of  the  surface  at  high  tempera- 
tures is  more  important.  Therefore  thermal  radiance 
of  the  surface  is  considered  in  the  analysis  of  the 
heat  balance  on  the  model  surface. 

Let  us  consider  a test  model,  which  is  exposed  to  a 
high  enthalpy  flow  field  (Fig.  3). 

The  heat  balance  on  the  model  surface  can  be  ex- 
pressed as 


Figure  3.  Heat  transfer  balance  at  the  model 
wall 

At  low  temperatures  of  the  model  in  a high  en- 
thalpy flow  field  the  thermal  radiance  (qr ) is  negli- 
gible compared  to  the  gas  heating  (Fig.  3).  Since 
the  difference  of  the  gas  and  surface  temperatures  is 
a driving  factor  of  the  heat  transfer  from  the  gas 
surface  to  the  model  surface  (eqs.  (1),  (2)  and  (3)), 
it  is  high  at  cold  surfaces  in  a high  temperature  gas 
flow.  Therefore  the  heat  conduction  into  the  model, 
i.e.  temperature  increase  in  the  model  structure,  is 
high.  Since  the  gas  heating  to  the  surface  and  ra- 
diation cooling  of  the  model  have  an  opposite  de- 
velopment with  increasing  temperature,  a so-called 
radiative  equilibrium  establishes  at  high  surface 


9A-6 


temperatures.  At  that  point  the  gas  heating  rate  and 
radiation  cooling  rate  are  equal.  The  heat  flux  rate 
to  a surface  in  radiative  equilibrium  can  be  deter- 
mined by  measuring  the  surface  temperature  and 
applying  the  equation  (11).  The  surface  temperature 
can  be  measured  with  non-intrusive  measurement 
techniques  like  pyrometer  or  IR-camera.  Besides  an 
accurate  temperature  measurement  a correct  value 
of  the  surface  emissivity  £ is  essential  for  the  de- 
termination of  the  so-called  hot  wall  heat  flux  rate. 
This  technique  can  only  be  applied  in  arc  heated  or 
induction  heated  facilities  with  long  testing  time 
(several  minutes),  where  radiative  equilibrium  can 
be  achieved. 


developments  are  the  boundary  conditions  of  this 
equation. 

Copper  slug  Copper  cylinder  Steel  cylinder 


In  order  to  determine  the  heat  flux  rate  to  models 
without  radiative  equilibrium  conditions  which  are 
achieved  in  long  duration  high  enthalpy  facilities, 
shock  tunnels  or  hot  shot  facilities  use  different 
measurement  techniques. 

4.1  Heat  storage  and  conduction  sensors 

4.1.1  Transient  heat  flux  probe  of  DLR 

The  transient  heat  flux  probe  designed  by  DLR  [4] 
consists  of  a copper  cylinder  with  two  embedded 
thermocouples  at  both  ends  of  the  cylinder  (Fig.  4). 
It  is  integrated  in  a water  cooled  probe  holder  al- 
most adiabatically. 

The  probe  is  swept  through  the  flow  field  at  con- 
stant flow  conditions.  Since  the  probe  temperature 
is  mostly  lower  than  500  K and  copper  surface  is 
catalytic,  this  probe  provides  cold  wall  heat  flux 
rate  to  a catalytic  surface.  The  heat  flux  rate  is  de- 
termined from  measured  temperature  development 
at  two  locations  and  applying  one  dimensional  heat 
conduction  equation.  The  measured  temperature 


Figure  4.  Transient  heat  flux  probe  of  DLR 

Based  on  the  description  in  the  Fig.  5 the  heat  con- 
duction equation  can  be  written  as  below: 


dL__K_  &T_ 
8t  pCp  dx 2 


(12) 


where  T,  t and  x are  the  temperature,  the  time  and 
the  axial  distance  from  the  front  of  probe  surface. 
The  parameters  K , p and  Cp  represents  the  thermal 

conductivity,  the  density  and  the  specific  heat  of  the 
copper  cylinder,  respectively. 


Ti 


(0.0-F, 


thermocouples 


Figure  5.  Mathematical  model  of  the  transient 
heat  flux  probe 
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The  initial  condition  is 

T(x,0)  = T0.  (13) 

The  boundary  conditions  are  given  by  two  func- 
tions F]  (?)  and  F2  ( t ): 

r(0,f)  = F1(t)  (14) 

T(L,t)  = F2(t ) (15) 

where  L is  the  length  of  the  cylinder. 

Since  the  solution  of  the  equation  (12)  requires 
complicated  Laplace  transformations,  it  is  not  de- 
scribed here  in  details.  Using  the  determined  tem- 
perature development  along  the  cylinder  by  solving 
equation  (12),  the  heat  flux  rate  to  the  front  surface 
of  the  probe  can  be  calculated  from  equation: 

(,6) 

Transient  heat  flux  measurement  technique  assumes 
that  no  heat  is  lost  by  conduction  in  radial  direction 
and  thermocouple  wires.  The  first  assumption  is 
generally  justified  by  keeping  the  contact  surface  of 
the  probe  to  the  probe  holder  as  small  as  possible 
(three  points  around  the  circumference  of  the  cylin- 
der and  use  of  an  insulation  material  as  interface 
surface).  Since  the  heat  conductivity  of  copper  is 
nearly  four  times  higher  than  that  of  Nickel  or 
Chromium  (NiCr-Ni  thermocouples  are  used)  [5] 
and  the  cross  section  area  of  thermocouples  is  very 
small  compared  the  cylinder  surface  area  the  heat 
sink  effect  of  thermocouples  is  small.  The  main  un- 
certainty of  this  technique  is  correlated  to  the  quali- 
ty of  vacuum  brazing  of  thermocouples  into  the 
copper  cylinder.  Even  small  cavities  at  the  meas- 


urement spot  of  thermocouples  in  the  cylinder  can 
influence  the  accuracy  of  the  probe  remarkably. 
Therefore  all  braze  spots  have  to  be  analysed  using 
X-ray  before  the  calibration  and  use  of  the  probe. 
The  temperature  dependency  of  the  material  prop- 
erties (p,  Cp,K ) is  considered  to  achieve  a good 

accuracy  of  the  measurement.  This  technique  also 
allows  to  measure  the  spatial  distribution  of  the 
heat  flux  rate. 

4.1.2  Heat  flux  sensor  of  NAL 

National  Aerospace  Laboratory  (NAL)  in  Tokyo 
has  developed  a heat  flux  sensor  for  high  tempera- 
ture application  using  tile  material  [6],  This  sensor 
was  successfully  applied  during  the  flight  of  the 
HYFLEX  vehicle  in  1996  [7].  It  consists  of  tile 
material  with  a high  temperature  coating  of  0.3  mm 
thickness  (Fig.  6).  The  sensor  is  attached  to  the  sub- 
structure of  the  vehicle  via  an  aluminium  flange 
and  screws. 


Figure  6.  NAL  heat  flux  sensor  [6] 
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Three  thermocouples  are  integrated  in  the  sensor  to 
measure  the  temperature  history  at  three  locations 
along  the  axis  of  the  plug  (Fig.  7).  The  first  ther- 
mocouple is  attached  to  the  rear  surface  of  the  high 
temperature  glass  coating  layer.  The  measured  tem- 
perature is  assumed  to  be  the  surface  temperature  of 
the  sensor.  The  much  higher  thermal  conductivity 
of  the  glass  coating  compared  to  the  ceramic  tile 
justifies  this  assumption.  Using  measured  tempera- 
tures with  three  thermocouples  as  boundary  conditi- 
ons under  the  assumption  of  two  dimensional  plane 
symmetry  and  applying  two  dimensional  heat  con- 
duction equation  the  heat  flux  rate  is  calculated. 

NAL  sensor  is  one  of  the  few  sensors,  which  can  be 
used  at  high  temperatures  up  to  1700  K for  flight 
experiments.  It  has  a response  time  of  about  one  se- 
cond. Since  it  is  made  of  tile  materials  of  the  vehi- 
cle, the  transferability  of  the  data  to  the  thermal  be- 
haviour of  other  TPS  components  of  the  vehicle  is 
more  reliable  compared  to  sensors  consisting  of 
metallic  materials.  Even  the  complex  data  reduction 
of  this  sensor  does  not  influence  its  advantages  for 
flight  application. 


Figure  7.  Principle  sketch  of  NAL  heat  flux 
sensor [7] 


4.2  Sensors  based  on  thick  wall  technique 

The  principle  of  thick  wall  heat  flux  sensors  is 
based  on  the  measurement  of  the  surface  tempera- 
ture of  a model  structure,  which  can  be  considered 
to  be  a semi-infinite  solid  (Fig.  8).  It  means  that 
during  the  measurement  time  period  the  conductive 
heating  does  not  reach  the  rear  surface  of  the  mo- 
del, i.e.  the  surface  temperature  at  that  location  re- 
mains constant. 


(Ut)=T0 


Figure  8.  Boundary  temperatures  of  a semi- 
infinite  model 


The  boundary  conditions  for  the  one-dimensional 
heat  conduction  equation  (eq.  (12))  would  be: 

r(0,f)  = F1(f)  = 7'(f)  (17) 


Qs(t)  = K 


(^(0) 

= K\ 

fdT(t)') 

V 5x  j 

x»0 

l dx  J 

(17’) 


T (L,t)  = T(co,t)  = T0 . 


(18) 


The  solution  of  the  equation  (12)  under  these 
boundary  conditions  leads  to  the  following  expres- 
sion for  the  heat  flux  rate  to  a semi-infinite  probe: 


q(t)  = 


pV 


T(t ) 1 r T(t)-T(x) 


- + 


7T 


VF  2 J (t  — x) 


,3/2 


dx 


•(19) 
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In  the  case  of  a constant  heat  flux  rate  equation  (19) 
takes  a simpler  form: 

<j=; j-^r  pw)-u  (20) 

This  technique  requires  a short  exposure  time  of  the 
sensor  to  the  flow.  The  duration  is  defined  by  the 
constancy  of  the  rear  surface  temperature  of  the 
slug.  Therefore  it  is  mainly  used  in  shock  tunnels, 
impulse  tunnels  or  hot  shot  facilities.  The  transition 
phase  of  the  flow  after  the  start  makes  the  heat  flux 
evaluation  in  these  facilities  difficult. 

4.2.1  Coaxial  surface  thermocouple 

The  principal  of  the  application  of  the  coaxial  sur- 
face thermocouple  for  the  heat  flux  rate  measure- 
ment is  the  measurement  of  the  surface  temperature 
of  a body,  which  can  be  considered  as  a semi-infi- 
nite solid.  One  thermoelectric  material  is  installed 
in  a second  tube  of  thermoelectric  material  concen- 
trically with  an  electrical  insulation  of  ceramic  base 
(Fig.  9).  The  hot  junction  of  the  thermocouple  is 
created  by  abrading  the  center  conductor  and  outer 
tube  together.  The  coaxial  thermocouple  assembly 
is  completed  by  attaching  thermocouple  lead  wires 
to  the  coaxial  thermoelements. 

The  heat  flux  rate  is  determined  from  the  measured 
surface  temperature  history  by  solving  equation 
(19)  or  (20)  using  a numerical  integration  techni- 
que. 

In  order  to  achieve  a representative  heat  flux  rate  to 
the  model  surface,  the  material  properties  of  the 

thermocouple  (i.e.  ^jpCpK  value)  has  to  match 
with  the  model  parameters.  This  technique  is  ap- 


plied in  short  duration  wind  tunnels,  continuous 
wind  tunnels  and  arc  jet  facilities.  The  main  ad- 
vantages of  the  coaxial  thermocouple  technique  are 
its  fast  response  time  (~  50  ps)  and  very  good  dura- 
bility. Since  the  output  is  self-generating,  no  cali- 
bration is  required.  It’s  application  is  restricted 
more  or  less  to  metal  models.  Weak  output  signal 
and  complex  data  reduction  are  the  disadvantages 
of  this  technique. 


COAX 

INSTALLED  IY  MESS  FITTING 

Mod  imUUTKHt 


Figure  9.  Sketch  of  coaxial  surface  thermo- 
couple [8] 

Based  on  already  available  coaxial  thermocouples, 
a coaxial  thermocouple  with  integrated  pressure 
gauge  has  been  developed  at  the  Technical  Univer- 
sity of  Aachen  (RWTH  Aachen)  (Fig.  10)  [9],  It 
allows  to  measure  the  temperature,  i.e.  heat  flux 
rate,  and  pressure  almost  at  the  same  location  (spa- 
tial deviation  between  two  sensors  is  about  1.5 
mm).  The  sensor  properties  are  matched  for  its  ap- 
plication in  short  duration  facilities. 
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Figure  10.  Coaxial  surface  thermocouple  with 
integrated  pressure  gauge  [9] 

4.2.2  Thin-film  gauge 

Thin  film  gauges  have  been  extensively  used  in  the 
aerodynamic  research  for  decades.  Mainly  a plati- 
num film  is  applied  in  a thickness  of  several  nano- 
meters on  the  substrate  by  sputtering,  vacuum 
deposition  or  painting.  Pyrex,  fused  Quarts  or 
MACOR  are  mainly  used  as  substrate.  Because  of 
its  high  positive  temperature-resistance  coefficient 
Platinum  provides  an  excellent  output  signal.  The 
temperature  data  deduced  from  the  measured  resis- 
tance change  is  reduced  to  the  heat  flux  rate  using 
semi-infinite  heat  flux  equation  (eq.  (19)  or  (20)). 
Figure  11  shows  a thin  film  gauge  developed  at 
RWTH  Aachen  [10].  Thin  film  gauges  are  very 
convenient  for  aerothermodynamic  application, 
since  they  have  a negligible  heat  capacity  compared 
to  the  substrate.  In  contrast  to  coaxial  surface  ther- 
mocouples thin  film  gauges  find  a broad  application 
on  non-metallic  materials.  The  other  remarkable 
advantage  is  their  easy  use  in  surface  curvatures  of 
small  radius.  The  response  time  is  very  short  (less 
than  5 ps).  But  these  sensors  are  very  sensitive  even 
to  small  particles.  The  calibration  is  time 


consuming.  The  influence  of  the  heat  conduction  on 
the  results  has  to  be  considered  in  the  evaluation  of 
the  data. 


Figure  11.  Thin  film  gauge  of  RWTH  Aachen 
[10] 

4.2.3  Null-point  calorimeter 

The  null-point  calorimeter  was  developed  by 
AEDC  [8].  Figure  12  shows  a sketch  of  this  sensor. 
A thermal  mass  with  a length  of  several  centimetres 
is  drilled  from  the  backside  to  a few  millimetre  dis- 
tance from  the  front  surface. 


[8] 


The  location  0,  b on  the  radial  centerline  of  the 
cylinder  cavity  is  defined  as  the  null  point.  The 
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measured  temperature  development  of  this  point  is 
assumed  to  be  identical  to  the  surface  temperature 
history  on  the  outside  surface  of  the  same  thermal 
mass  without  cavity.  In  order  to  justify  this  as- 
sumption the  geometry  of  the  sensor  is  specified  by 
a very  detailed  thermal  analysis.  According  to  null 
point,  analysis  states  that  the  temperature  history  at 
the  null  point  would  be  very  close  to  the  surface 
temperature  history  in  the  absence  of  the  hole,  if  the 
ratio  of  the  hole  radius  to  the  axial  distance  of  the 
point  about  1.4.  Therefore,  the  measured  null  point 
temperature  with  a Chromel-Alumel  thermocouple 
can  be  inserted  into  the  equation  (19)  to  determine 
the  heat  flux  rate  to  a semi-infinite  mass.  The  ther- 
mocouple wires  are  normally  attached  by  vacuum 
brazing.  As  it  is  shown  in  Fig.  13  the  null-point 
body  has  a slight  chamber  at  the  top  and  bottom, 
which  creates  an  effective  circumferential  dead  air 
space  along  the  length  of  the  cylinder  to  enhance 
one-dimensional  heat  conduction  and  prevent  radial 
heat  transfer. 


Figure  13.  Cross  section  view  of  null-point  ca- 
lorimeter [8] 

This  sensor  is  mainly  used  at  the  stagnation  point 
test  configuration  of  arc  heated  facilities.  Its  re- 
sponse time  is  1-2  ms.  Very  high  heat  flux  rates  up 


to  284  MW/m2  can  be  measured  with  this  sensor. 
Shortcomings  of  the  null-point  calorimeter  are  the 
difficulty  in  the  achievement  of  consistency  in 
thermocouple  wire  attachment,  difficult  calibration 
procedure  and  its  limitation  to  high  heat  flux  rates 
at  short  time  duration. 

4.3  Calorimetric  sensors 

4.3.1  Thin-skin  technique 

The  operation  principle  of  the  thin-skin  technique 
involves  measuring  the  slope  of  the  back  surface 
temperature  history  of  a thin  skin  model. 

The  heat  flux  rate  can  be  calculated  using  the  rela- 
tion: 

<’=pc's^  (2I) 

where  s is  the  thickness  of  the  thin  wall. 

This  technique  assumes  that  the  heat  conduction  to 
adjacent  structures  is  negligible  during  the  meas- 
uring time.  The  exposure  time  of  the  sensor  has  to 
be  specified  in  dependence  on  the  material  proper- 
ties and  heat  flux  rate  level,  i.e.  temperature  gradi- 
ent on  the  back  surface.  Simple  data  reduction  is  a 
big  advantage  of  this  technique.  But  the  influence 
of  the  lateral  heat  conduction  and  temperature  de- 
pendency of  the  material  properties  on  the  meas- 
urement accuracy  is  a problem.  Therefore  the  fabri- 
cation of  a thin-skin  model  is  expensive. 

4.3.2  Slug  calorimeter 

The  slug  calorimeter  technique  is  based  on  the  stor- 
age of  all  the  heat  which  enters  upon  the  sensor 
surface  in  a calorimetric  mass  during  the  measuring 
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time.  A metallic  disc  is  mostly  used  as  calorimetric 
mass.  The  slug  is  integrated  in  a thermal  insulator 
to  avoid  radial  heat  losses  (Fig.  14). 


calorimetric  mass 


Figure  14.  Sketch  of  slug  calorimeter 

Again  equation  (21)  is  applied  for  the  calculation  of 
the  heat  flux  rate  from  the  measured  temperature 
history  on  the  backside  of  the  disc.  The  accuracy 
depends  strongly  on  the  quality  of  the  insulation, 
i.e.  prevention  of  radial  heat  conduction. 


The  water  mass  flow  rate  and  the  temperature 
increase  of  the  water  are  measured  precisely.  The 
heat  flux  rate  is  calculated  using  following  relation: 


^cCpc(7-2c-7-1c) 

A 


where 


(22) 


A : cross  section  area  of  the  probe 

Cpc : specific  heat  of  water 

mc : water  mass  flow  rate 

Tlc : water  inlet  temperature 

T2c  : water  outlet  temperature 


Material  of  different  catalysis  can  be  used  as  probe 
material.  DLR  and  VKI  have  manufactured  calori- 
meter probes  based  on  IPM-concept  (Fig.  16). 


4.3.3  Water  cooled  calorimeter 

The  Institute  of  Mechanical  Problems  (IPM)  in 
Moscow  developed  a water  cooled  calorimeter  to 
measure  the  heat  flux  rate  in  steady  state  flow  con- 
ditions [11].  This  technique  is  based  on  the  meas- 
urement of  the  energy  absorbed  by  cooling  water, 
which  removes  the  incident  heat  from  the  back 
surface  of  the  thin  metallic  disc  of  a hollow  cylin- 
der (Fig.  15). 


removable  cylinder 


Figure  16.  Water  cooled  calorimeter  probe  of 
DLR  based  on  EPM  concept 


Figure  15.  Sketch  of  water  cooled  steady  state 
calorimeter 


4.4  Heat  conduction  gauges 

4.4.1  Gardon  gauge  (asymptotic  calorimeter) 

Gardon  gauges  consist  basically  of  a thin  metallic 
foil  suspended  over  a cavity  in  a heat  sink.  It  is 
thermally  and  electrically  attached  to  the  heat  sink 
at  the  periphery  of  the  cavity  using  metallurgical 
bonding  techniques  (Fig.  17).  The  metallic  foil  as 
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the  first  thermoelectric  material  (mostly  constantan) 
creates  together  with  to  the  second  thermoelectric 
material  of  periphery  (copper)  one  junction  of  the 
thermocouple.  A wire  of  second  thermoelectric  ma- 
terial is  attached  to  the  rear  surface  center  of  the 
foil  to  form  a differential  thermocouple.  It  measures 
the  temperature  difference  between  the  center  and 
the  edge  of  the  foil. 

The  one-dimensional  heat  conduction  in  radial  di- 
rection can  be  described  in  polar  co-ordinates  as: 

pcp  dT  1 8T  82T  q 

= 1 — h (23) 

K dt  r dr  dr 2 s./C 

with  boundary  conditions 

T (r,0)  = Ts  at  0 <r<R  (24) 

and 

T (R,  t)  = Ts  for  0<f<oo.  (25) 

The  steady  state  solution  of  the  equation  (23)  under 
these  boundary  conditions  is: 

A<-ic 

<26) 

The  heat  flux  rate  at  the  foil  center  ( r = 0 ) would 
be 

<?  = ^(r-n).  (27) 

Equations  (26)  and  (27)  show  that  the  temperature 
difference  between  the  center  and  edge  of  the  foil  is 
directly  proportional  to  the  heating  rate  q . A fur- 
ther thermocouple  is  usually  integrated  inside  the 


cylinder  to  provide  absolute  gauge  temperature  for 
the  calculation  of  heat  transfer  coefficient  in  aero- 
dynamic application.  The  influence  of  the  heat  sink 
temperature  should  be  considered  by  the  specifica- 
tion of  the  sensor.  Copper-Constantan  combination 
meets  the  requirements  with  respect  to  the  linearity 
of  the  sensor. 


constantan  foil 


Figure  17.  Principle  scheme  of  Gardon  gauge 

Gardon  gauges  have  a broad  application  in  wind 
tunnel  heat  flux  measurements  and  heat  transfer 
determination  during  flight.  Uncooled  sensors  can 
only  be  applied  at  surface  temperatures  up  to 
1 100°C.  The  response  time  of  these  sensors  is  about 
1 second.  Since  the  heat  flux  rate  is  directly  pro- 
portional to  the  temperature  difference,  the  data 
reduction  is  very  easy.  Low  sensitivity  and  limited 
surface  temperature  are  main  disadvantages.  Gar- 
don gauge  signal  in  high  enthalpy  facilities  is  noisy. 

4.4.2  Thin  thermopile  sensors 

The  thermopile  measurement  technique  is  based  on 
the  measurement  of  the  temperature  change  across 
the  thermal  resistance  elements,  which  are  depos- 
ited on  a heat  sink  material.  The  thermopile  consists 
of  multi  thermocouple  pairs  to  achieve  measurable 
voltage  drop  across  the  sensor. 
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Vatell  Corporation  has  developed  a micro-sensor 
(HFM)  based  on  this  technique.  Two  different  thin 
films  are  deposited  on  the  aluminium  nitride  sub- 
strate of  HFM.  A thin  film  of  less  than  2 microns  of 
the  heat  flux  sensor  HFS  allows  to  achieve  very  fast 
response  of  about  10  microseconds.  The  heat  flux 
sensor  is  surrounded  by  a temperature  sensor  made 
of  platinum  resistors  (RTS).  Temperature  is  deter- 
mined by  passing  a small  constant  current  through 
the  resistance  and  measuring  the  resulting  voltage. 
Substrate  surface  temperature  is  used  to  consider 
the  change  of  the  material  properties  with  the  tem- 
perature in  the  heat  flux  determination.  It  may  also 
be  used  to  check  the  calibration  of  the  HFM  sensor 
and  measure  the  heat  transfer  coefficient  in  aerody- 
namic application. 

In  order  to  use  the  linear  range  of  the  temperature 
dependency  of  substrate  material  properties,  the 
heat  transfer  to  the  surrounding  model  is  forced  by 
using  nickel  or  copper  housing,  which  is  bonded  to 
the  substrate  with  high  grip.  The  surface  of  the  thin 
films  can  be  coated  with  layers  of  different  emis- 
sivity,  which  is  very  essential  for  the  calibration  of 
the  sensor  using  radiation  sensors.  HFM  sensors 
can  be  used  at  temperatures  up  to  1 125  K without 
cooling. 

In  the  frame  of  X-38  programme  OHB  will  use  a 
HFM  sensor  integrated  in  a flexible  insulation  ma- 
terial (FEI)  and  perform  heat  flux  measurement  on 
the  Lee  side  of  the  vehicle  during  the  first  flight  of 
the  X-38  demonstrator  (Fig.  18)  [12]. 

The  OHB  set-up  including  HFM  sensor  has  been 
qualified  at  flight  conditions  of  X-38  demonstrator 
in  the  arc  heated  facility  LBK  of  DLR  (Fig.  19) 
[13]. 


Ceramic  Plug 


Figure  18.  HFM  sensor  integrated  in  the  OHB 
flight  component  [12] 


Figure  19.  Qualification  test  on  OHB  flight 
component  with  a HFM  sensor  in 
LBK 


5.  Calibration  of  heat  flux  sensors 

Calibration  of  heat  flux  sensors  is  a sophisticated 
work  for  engineers.  Since  heat  sources  based  on 
convective  heating  is  difficult  to  handle  with  re- 
spect to  their  repeatability,  mostly  radiative  cali- 
bration sources  are  used.  AEDC  has  developed  two 
such  calibration  sources.  The  first  source  uses  a 1 .6 
kW  xenon  arc  lamp  focussed  on  a small  area  to 
achieve  high  heat  flux  rates  up  to  28.4  MW/m2.  To 
improve  the  uniformity  of  the  radiant  intensity, 
which  is  a general  problem  of  most  of  the  radiant 
heat  sources,  an  integrator  consisting  of  an  alu- 
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minium  tube  with  a highly  polished  inside  surface 
was  installed  with  its  entrance  at  the  focal  plane  of 
the  calibrator.  Another  portable  calibrator  for  lower 
heat  flux  rate  up  to  1 .4  MW/m2  has  a 200  W tung- 
sten filament  lamp  as  heat  source.  The  radiant  en- 
ergy of  the  lamp  is  focussed  on  a small  spot  of 
about  7.6  mm  at  one  of  the  focal  point  of  the  ellip- 
soidal reflector  with  polished  surface  (Fig.  20).  The 
sensor  is  placed  at  this  focal  point. 

Both  calibration  techniques  require  a reference 
sensor  in  addition  to  the  sensor,  which  has  to  be 
calibrated. 

The  calibration  unit  of  Vatell  uses  the  radiation  of  a 
electrically  heated  thin  graphite  plate  as  heat  source 
[12].  The  reference  and  calibration  sensor  are  in- 
stalled at  each  sides  of  the  plate  and  heated  with 
same  radiance  intensity.  In  order  to  avoid  any  con- 
vection and  oxidation  processes,  the  complete  fur- 
nace is  set  evacuated  down  to  1 hPa  pressure  level. 
Heat  flux  rate  up  to  3.4  MW/m2  can  be  achieved  in 
this  heat  flux  calibration  source. 


6.  Comparative  heat  flux  measurements  in  LBK 

Transient  heat  flux  probe  and  water  cooled  steady 
state  calorimeter  of  LBK  (Figs.  4 and  16)  are  used 
for  the  measurement  of  cold  wall  heat  flux  rate  to  a 
copper  surface  at  same  test  conditions.  Tests  were 
performed  in  the  stagnation  point  test  configuration 
of  LBK.  Both  flow  conditions  were  chosen  in  such 
a way  that  moderate  heat  flux  rates  can  be  set.  A 
SiC  model  with  the  same  geometry  as  heat  flux 
probes  was  also  tested  to  determine  the  hot  wall 
heat  flux  rate.  The  duration  of  the  tests  on  the  SiC 
probe  was  long  enough  to  reach  radiative  equilib- 
rium. As  mentioned  before  the  hot  wall  heat  flux  is 
determined  from  the  measured  surface  temperature 
with  a pyrometer.  Measured  cold  and  hot  wall  heat 
flux  rates  are  listed  in  Table  1. 


Test  condition 

I 

II 

qCft  [kW/m2] 

725 

1080 

qc„  [kW/m2] 

680 

1210 

qh  [kW/m2] 

303 

351 

SWT  HEATED  HOLDER 


SfOT  HEATED 


SCHMIDT-WELTED 
GA6E  STAND/UtD 

SLOT  FOR  SMUTTED 


- INTERCHANGEAIIE  SENSOR/MODEL 
HOLDER 

10*  HAlf  ANGLE  SRHERE 
MODEL 


NULL-POINT  CALORIMETER 
SHUTTER 

INTERCHANGEAIIE 
STANDARO/SENSOI 
HOLDER 


Table  1.  Cold  wall  heat  flux  rates  to  a copper 
surface  measured  using  transient 
and  calorimeter  probes  in  compari- 
son with  hot  wall  heat  flux  rates  to  a 
SiC  model  in  LBK 


Although  at  test  condition  I the  agreement  between 
cold  wall  heat  flux  rates  measured  with  the  tran- 
sient and  calorimeter  probes  ( qc.i  and  qc.s , respec- 


Figure  20.  Sketch  of  portable  calibrator  of 
AEDC  [8] 


tively)  is  acceptable,  a deviation  of  around  1 1 % is 
noticed  at  the  test  condition  II.  A systematic  study 
including  a detailed  parameter  variation  is  under- 
way to  clarify  this  discrepancy.  Another  very  im- 
portant results,  which  is  strongly  related  to  non- 
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equilibrium  effects  of  high  enthalpy  flow  on  model 
surfaces  with  different  catalycity,  is  the  difference 
between  cold  and  hot  wall  heat  flux  rates.  As  it  can 
be  seen  in  the  Table  1 the  cold  wall  heat  flux  rate  is 
more  than  twice  of  the  hot  one.  Since  the  gas  tem- 
perature at  the  stagnation  point  in  the  boundary 
layer  flow  of  the  model  of  more  than  4000  K is 
remarkably  higher  than  the  model  surface  temper- 
ature and  the  difference  in  the  surface  temperature 
of  SiC  and  copper  cylinder  of  cold  wall  heat  flux 
probes  is  about  1500  K,  this  big  difference  has  to  be 
caused  by  additional  gas  surface  interaction  effects. 
At  temperatures  beyond  2500  K molecular  oxygen 
in  the  flow  field  of  LBK  is  completely  dissociated. 
Spectroscopic  measurements  and  numerical  results 
state  that  the  flow  field  in  the  free  stream  and 
behind  bow  shock  ahead  of  the  model  is  in  thermo- 
chemical non-equilibrium  (nearly  frozen).  There- 
fore the  surface  catalysis  has  an  significant  influ- 
ence on  the  heat  flux  rate.  It  is  well  known  that  at 
above  specified  test  condition  a SiC>2  layer  develops 
due  to  passive  oxidation  on  the  SiC  surface.  This 
leads  to  a low  catalytic  surface.  In  contrary  copper 
surface  is  nearly  fully  catalytic.  Therefore  strong 
exothermic  recombination  reaction  takes  place  on 
the  copper  surface  and  enhance  the  heat  flux  rate. 
Especially  non-equilibrium  effects  on  the  gas  sur- 
face interaction  keeps  heat  flux  measurements  in 
high  enthalpy  flow  fields  very  challenging  and  re- 
quire further  improvements. 

7.  Concluding  remarks 

Parallel  to  the  progress  in  material  science,  elec- 
tronics, and  spectroscopic  techniques  the  accuracy 
and  application  range  of  heat  flux  sensors  has  been 
improved  in  the  last  decades.  But  compared  to  other 


parameters  in  aerothermodynamics  like  pressure, 
flow  rate,  etc.,  where  accuracy  and  repeatability 
values  are  specified  within  1%,  usually  an  accuracy 
value  within  5%  cannot  be  guaranteed  for  heat  flux 
sensors,  yet.  Another  challenge  for  engineers  is  the 
development  of  sensors  of  high  accuracy  for  high 
temperature  application  in  a reactive  environment. 
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